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In recent years the precision of lattice calculations has improved hugely, and the results are mak- 
ing a very significant impact in particle physics phenomenology. Indeed there is no alternative 
general method which can be used in the evaluation of nonperturbative strong interaction effects 
for a wide variety of physical processes. In this talk I discuss a selection of topics in flavour 
physics, including mature quantities for which lattice calculations have been performed for a long 
time (e.g. the determination of the V us CKM matrix element and Bk), quantities which we are 
now learning to study (e.g. K —> nn decays amplitudes and the spectrum and mixing of 77 — rj' 
mesons) and important phenomenological quantities for which a large amount of experimental 
data is available but which we do not yet understand how to approach in lattice simulations (e.g. 
nonleptonic B-decays). The improvement in precision and the extension of the range of processes 
which can be studied using lattice QCD has to be continued vigorously if precision flavour physics 
is to play a complementary role to large P± discovery experiments at the LHC in unravelling the 
next level of fundamental physics. 
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1. Introductory Comments 

Recent progress in lattice phenomenology has been truly impressive, due to improved algo- 
rithms and theoretical developments as well as more powerful computing resources. It is clearly 
impossible to do justice to the title of this talk by discussing all the exciting new developments and 
results and I will not try to do this. At this conference there have been a number of excellent review 
talks on specific topics including those in refs. [0-01 as well as many parallel sessions in which 
the new ideas and results have been presented explicitly. I will therefore not attempt to make a 
systematic compilation of the latest lattice results. I will also not discuss the different formulations 
of lattice QCD or effective theories of heavy quarks. Instead I will discuss and give examples of 

1. Mature quantities, i.e. ones which have been studied successfully and with improving preci- 
sion for many years. These include the determination of the CKM matrix element V us and the 
Bk parameter of neutral kaon mixing, which I will discuss. The natural question which arises 
for such quantities is what next? Does the phenomenology require us to strive towards even 
higher precision or should we set our priorities towards the evaluation of other quantities? 

2. Continuing attempts to extend the range of physical quantities for which lattice simulations 
can contribute to the quantitative control of the nonperturbative QCD effects. In this context 
I will discuss attempts to evaluate K — » %% decay amplitudes and a recent study of the 17 — 17' 
system. I will also briefly mention relatively early attempts to evaluate long-distance effects, 
encoded in the matrix elements of time-ordered (non-local) products of operators. 

3. Quantities for which we don't yet know how to formulate the lattice calculation of the corre- 
sponding non-perturbative QCD effects. For illustration I will discuss two-body nonleptonic 
B-decays, a hugely important set of processes for which there is a large amount of experimen- 
tal data and yet lattice simulations are not playing any significant role in the phenomenology. 

The examples I use in this talk are taken from Flavour Physics, but of course we also need to re- 
member the contributions which lattice QCD is making to the development of our understanding 
of hadronic structure and to the determination of quark masses and the coupling constant a s . The 
mission of precision flavour physics is to play a complementary role to large p± experiments in dis- 
covering and unraveling the next layer of fundamental physics beyond the standard model. If, as is 
expected, or at least hoped, the LHC experiments discover new elementary particles then precision 
flavour physics will be necessary to determine the underlying theoretical framework and lattice 
simulations will be central to this endeavour, with the specific task of quantifying nonperturbative 
hadronic effects. The discovery potential of precision flavour physics also should not be under- 
estimated however, with a real possibility that current tensions within the standard CKM analysis 
will be confirmed and become inconsistencies which will have to be explained by new physics. 
Even restricting the discussion to flavour physics, I have to be selective and have naturally chosen 
to discuss the status of some topics which have been of direct interest to me in recent years, and 
which I believe are indicative of the exciting progress and prospects for the subject more widely. 

Before proceeding to the discussion of specific processes, I start with some comments about 
the use of continuum perturbation theory which is necessary to relate renormalization schemes 
which can be simulated in lattice calculations and ones based on dimensional regularization (such 
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as MS) which cannot, but which are generally used in the perturbative evaluation of Wilson co- 
efficient functions. The main point of this discussion is to underline that the precision of current 
lattice computations is such that, in order to get the maximum scientific benefit from the results 
we need to work with the higher-order perturbation theory community to obtain Wilson coefficient 
functions in schemes we can simulate. After this digression I proceed to a discussion of the status 
of the evaluation of a variety of physical quantities, starting with standard ones such as V us and Bk, 
through quantities which we are beginning to evaluate reliably such as K — > %% decay amplitudes 
and 17 — 77' masses and mixing and finally to some quantities which we currently do not know how 
to evaluate at all (non-leptonic B-decay amplitudes). 

1.1 Continuum Perturbation Theory 

Lattice simulations are used to compute the long-distance non-perturbative effects in quan- 
tities such as the matrix elements of local composite operators or the QCD parameters a s (the 
strong coupling constant) and the quark masses. These quantities require renormalization and so 
have to be combined with perturbative calculations before they can be used in predicting physical 
observables. This is sketched in the following oversimplified picture: 



where C (perhaps a Wilson coefficient function) contains the short-distance physics and is calcu- 
lated in perturbation theory and the long-distance effects are contained in the matrix element(s) 
of one or more local operators. The renormalization scheme and scale dependence cancels be- 
tween the two factors C and ( / 1 \ i) provided that they are calculated in the same renormalization 
scheme. Because of its practical advantages C is usually calculated in a scheme based on dimen- 
sional regularization, such as the MS scheme. Unfortunately however, we are not able to perform 
lattice simulations in a non-integer number of dimensions and hence MS, while being the scheme 
of choice in continuum perturbative calculations, is not directly useful in combining the perturba- 
tive and lattice results. In particular, it is not possible to calculate the matrix element in the MS 
scheme entirely non-perturbatively. 

One possibility for matching C and the matrix element, as long as the ultraviolet cut-off 
{aT x , where a is the lattice spacing) and the renormalization scale (/i) are both sufficiently large, is 
to use perturbation theory to relate (/ 1 Mt (a) | i), computed in the bare lattice theory and with bare 
operators latt (a), to the corresponding renormalized operators MS (/i). In addition to the techni- 
cal difficulty of performing perturbative calculations in the lattice theory beyond one-loop order, it 
is found however, that the series frequently converge poorly. It is therefore preferable to perform 
Non-Perturbative Renormalization (NPR), by choosing a renormalization condition which can be 
imposed directly in lattice computations [Q| and examples of such schemes include the momentum 
schemes RI-Mom [|J and its recent generalization RI-SMom [[7], ||] and schemes based on the use 
of the Schrodinger Functional Having obtained the matrix element in such an intermediate 
scheme, we need to combine it with C calculated in the same scheme or equivalently to translate 
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the matrix element from the intermediate scheme to MS. The precision of lattice calculations is 
now such that this translation from the intermediate scheme to MS needs to be performed beyond 
one-loop order and hence is left most effectively to the professional N n LO perturbation theory 
specialists. It is therefore pleasing to see such calculations being performed; this enhances con- 
siderably the phenomenological reach of the lattice calculations. In addition to John Gracey, who 



has been performing such calculations for some time (see for example |]10p), there have been recent 
calculations relating the quark masses in the RI-SMom and MS schemes at two-loop order [O, O]. 
Note also the joint work by the HPQCD lattice collaboration and the Karlsruhe perturbation theo- 
rists on the evaluation of the charm-quark mass by matching moments of correlators on the lattice 



and in the continuum [13] 



The principal lesson of this section is that close collaboration with the QCD perturbation the- 
ory community is increasingly necessary to optimise the significance of our lattice results. 

2. Vus from Lattice simulations 

Lattice calculations of the decay constants f n and fx are by now standard and from the ratio 
Ik/ fit, combined with the experimental leptonic widths, we obtain the ratio V us /V u d. The current 
status obtained from simulations with Nf = 2 + 1 flavours of sea quarks is [ fl4| ] 

^ = 1.193 ±0.006. (2.1) 

J7t 



In the last 6 years or so, following the suggestion of Becirevic et al. Jl 5fJ , it has also become 
possible to determine V us precisely by combining the experimental results for K — > % semileptonic 
decays with lattice determinations of the form factor / + (0) = f°(0), where the argument in 
the parentheses indicates that the four-momentum transfer q between the kaon and pion satisfies 
q 2 = [0- 18]. The FLAG collaboration summarises the current status of the results for the form 



factor as [14] 



/+(0) = 0.956 ±0.008. (2.2) 

For both the quantities fx/fn and / + (0), the calculational techniques are such that we would obtain 
precisely 1 in the SU(3) flavour symmetry limit (m u = m c j = m s ) and so it is the difference from 
one which we are actually computing. The Ademollo-Gatto theorem implies that these corrections 
are small for / + (0), and the main uncertainty for this quantity is due to the chiral extrapolation. 



The precision of the results in eqs. (2.1) and (2.2) is truly remarkable when compared to what was 
possible just a few years ago. 

The recent results are nicely summarized in figure [j], also from FLAG [14], in which are 
shown: 

(i) the allowed regions from calculations of fx/fn with Nf = 2+1 and Nf = 2 flavours of 
sea-quark (shown separately); 

(ii) the allowed regions from calculations of / + (0) with Nf = 2 + 1 and Nf = 2 flavours of 
sea-quark (shown separately); 

(iii) the allowed regions combining the calculations of fx/fn and / + (0) with Nf = 2 + 1 (red 
oval) and Nf = 2 (dashed blue oval). 
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Figure 1: The plot compares the information for \V uc i\ and |V" MS | obtained on the lattice with the experimental 
result extracted from nuclear j3 transitions. The dotted arc is part of the circle \V u d\ 2 + \V US \ 2 — 1, representing 
the correlation between \V llc i\ and \V US \ that follows if the three-flavour CKM-matrix is unitary, [pft. 

Also shown on this plot is the result for V lu / from super-allowed nuclear j3 -decays and an arc of 
the circle |V„d| 2 + \ V US \ = 1. Since \V u b\ 2 is very small compared to the uncertainties in \V u d\ 2 and 
\Vus\ 2 , this circle represents the unitarity condition on the first row of the CKM matrix. Although 
it would have been more exciting to find an inconsistency, we see that everything is remarkably 
consistent with the standard model, a point I underline further in the next subsection. 

2.1 V U s within the Standard Model 



I'd now like to share a simple analysis, which I learned from my colleagues in the FLAG [ |14fJ 
collaboration, which determines V us in the Standard Model without lattice calculations and which 
underlines again how remarkably consistent the lattice results for fx/ fx and / + (0) are with Stan- 
dard Model expectations. Experimental studies of leptonic decays of kaons and pions give us two 
precise constraints for the four quantities, V u d, V us , fx/fn and / + (0) Jl9|]: 

^^=0.27599(59) and |V^/ + (0)| =0.21661(47). (2.3) 

V u d J 71 

Within the standard model, a third equation is provided by the unitarity constraint: 

\V ud \ 2 + \V us \ 2 = l, (2.4) 
where |V„/,| 2 has been dropped from the left-hand side because it is smaller than the uncertainties 



in the remaining terms. Now eqs. (2.3) and (2.4) give 3 equations for four unknowns. As the fourth 



equation we can take the value of V uc i from the recent analysis in ref. [ |20| | based on 20 different 
superallowed nuclear j8 -transitions: 

\Vud \ =0.97425(22). (2.5) 
Combining equations ( |2.3| ), ( ^4fj and ( |2"3| ) we obtain: 

\Vus\ =0.22544(95), / + (0) = 0.9608(46), ^ = 1.1927(59), (2.6) 

Jn 



in excellent agreement with the lattice results for fx/ fn and / + (0) in (2.1) and (^2f) 
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m„ 


/°(<?Lx) 


670 MeV 


1.00029(6) 


575 MeV 


1.00016(6) 


555 MeV 


1.00192(34) 


470 MeV 


1.00272(34) 


415 MeV 


1.00887(89) 


435 MeV 


1.00416(43) 


330 MeV 


1.02143(132) 


375 MeV 


1.00961(123) 






300 MeV 


1.01923(121) 






260 MeV 


1.03097(224) 



Table 1: Values of f(q 2 max ) as a function of m n obtained by the RBC-UKQCD collaboration with 2+1 



flavours of Domain Wall Fermions (left two columns) [ 16 1 and by the ETM collaboration with 2 flavours of 
twisted mass fermions (right two columns) [|l7|]. 



2.2 K(3 decays and SU(2) chiral perturbation theory 

I end this section with some observations relating to semileptonic K — > 71 decays (Ka decays) 
using SU(2) chiral perturbation theory. Lattice calculations of / + (0) = /°(0) start with the eval- 
uation of /°(<7max)> where g^ ax = [mx — m K ) 2 is the largest physical value of q 2 (where q is the 
momentum transfer) and corresponds to the kinematics in which the initial kaon and final pion are 
both at rest. /°(^max) can De calculated with excellent precision as illustrated in table RL where 
the results from the RBC-UKQCD and ETM collaborations are presented as functions of the pion 
mass. The reason that I present these results here, is that the entries in table [T] appear to a long 
way below the value expected in the SU(2) chiral limit. This value is given by one of the very few 
analytic non-perturbative results in QCD, the Callan-Treiman relation, which I present here in the 
form 

/°(<7Lx) = 7^1-26 (2.7) 

Jit 

where all the quantities are evaluated in the SU(2) chiral limit, m u = raj = . It is a little puzzling 
that while the values of /°(<?max) m table [I] are increasing as m n decreases, they are only doing so 
very slowly. In ref. [pT]|, we used SU(2) chiral perturbation theory to investigate whether the very 
slow increase in /°(<7m ax ) observed in table [I] as m % decreases is expected to accelerate towards the 
value in eq. ( |2"77| ) for smaller values of the pion mass. We found that the one-loop chiral logarithms 
have a large coefficient and are of the correct size to account for the difference but they have the 
wrong sign, implying that the analytic terms (or higher order chiral logarithms) should account 
for approximately twice the difference between the results in table [l] and eq. (2.7). Of course the 



analytic terms (both linear and quadratic in m K ) are proportional to unknown low-energy constants 
(LECs) and hence are not calculable. In SU(3) chiral perturbation theory on the other hand, at 



one-loop order the LECs can be expressed in terms of fx/fx Q22| , [23j. Estimating the SU(2) LECs 
by converting results from SU(3) ChPT suggests that the analytic terms have the correct sign and 
(large) magnitude to account for the difference. 

It should be stressed here that since, apart from the values of the LECs, the above discussion 
relied only on SU(2) chiral perturbation theory and therefore the same features hold for B — > % and 
D — > 71 semileptonic decays. 
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Publication 






N f 


B K 


RBC-UKQCD 2007 
Aubin et al. 2009 [ 
Bae et al. 2010 [| 
RBC-UKQCD 2010 


[ 

2 
\ 
\ 


T§ 

8] 
] 

30] 


2+1 
2+1 
2+1 
2+1 


0.720(13)(37) 
0.724(8)(29) 
0.724(12)(43) 
0.749(7)(26) 


JLQCD 2008 [f31 
ETMC 2010 [32 


I 




2 
2 


0.758(6)(71) 
0.729(30) 



Table 2: A comparison of recent results for the renormalization group invariant Bk . Nf denotes the number 
of dynamical quark flavours. In each case, the first error is statistical and the second is systematic. 



Of course, we really want to know the chiral behaviour of the form-factor at the standard 
reference point q 2 = 0. At this point, the energy of the pion in the rest frame of the kaon is 
approximately mk/2 and so cannot be considered soft in SU(2) chiral perturbation theory. In spite 
of the hard external pion, we found that it is possible to evaluate the chiral logarithms, since they 



can be calculated from soft internal loops, finding y21| ] 

/°(0)=/+(0) 

r(o) 



3 


2 


log 


(ml 


4 


167T 2 / 2 


U 2 


3 


2 


log 


( m l 


4 


\6ll 2 f 2 


U 2 



+ c + m. 



(2.8) 
(2.9) 



where F± and c± are unknown LECs. In this way we have some information about the chiral 
behaviour of the form factors at q 2 = 0. Since the chiral extrapolation of the lattice results is 
the largest uncertainty in f° (0) any information about the chiral behaviour is useful (it would be 
very useful indeed if the results in eqs^l^) and could be extended to two loops in chiral 
perturbation theory). This hard-pion chiral perturbation theory approach has been generalized to 
K — > kk decays ||24|] and to B — > n and D — > n semileptonic decays [25]. 



3. Bk 

One of the very important quantities in particle physics phenomenology for which the pre- 
cision of lattice results has improved hugely in recent years has been Bk, the bag parameter of 
neutral kaon mixing. It is defined as the suitably normalized matrix element of the AS = 2 oper- 
ator (JlY^sl) (JlY^sl) between an initial K and final ^° state (where L denotes the left-handed 
component of the spinor field). As recently as five years ago, Chris Dawson, the rapporteur at 
this conference, was reporting results with a 12% error largely due to the fact that there were in- 
sufficient unquenched results at light masses A detailed review of recent results for Bk has 
been presented at this conference by Jack Laiho [Sj, for illustration here I tabulate the key results 
in table ^. Not only are the errors reduced by a factor of 3 or so , but the results in table || were 
obtained with a variety of fermion actions and techniques adding significantly to our confidence in 
the evaluation of the systematic uncertainties. 



7 



Lattice Phenomenology 



Christopher Sachrajda 



Collaboration(s) 


Re Ar/Re A 2 


e'/e 


RBC [ 
CP-PAC 
Experin 


41] 
lents 


25.3± 1.8 
9-^12 
22.2 


-(4.0±2.3) x 1(T 4 

(-7-h-2)xl0~ 4 
(17.2±1.8) x 1(T 4 



Table 3: Quenched 2001 results on the AI = 1/2 rule and e'/e obtained from K — > % and K — > vacuum 
matrix elements using the lowest order term in the SU(3) chiral expansion. 



Lattice results for Bk are an important ingredient in global studies of the unitarity triangle. 
Although in general the remarkable consistency of the information from different physical pro- 
cesses significantly restricts the possible parameter space for new physics, a number of tensions 
have arisen in recent years at the 1.5-3 standard deviation level (see for example [33- 37]), and 
the results for Bk contribute to these tensions. Calculations of Bk will continue towards 1% preci- 
sion, but already with current precision it is necessary to begin considering corrections which had 
previously been neglected. For example it has been stressed that in the theoretical expression for 
£k, the parameter monitoring indirect CP-violation in K — > Tin decays, we should include terms 
proportional to Im Ar/Re Ao (where Ao is the K — )■ Tin amplitude with the two pions in a state with 
isospin 0) and recognise that the phase arctan(2AM^/Ar) is not precisely equal tt/4 (AMk and Ar 
are the differences of the masses and widths of the Ki and K$ mesons) [ 3 8 - 4C ] . Necessary theoret- 
ical improvements include the evaluation of long distance effects, for which we need to determine 
matrix elements of the time-ordered product of two operators (see sec. ^). 



4. K-^txtx decays 

An important challenge for the lattice phenomenology community is the reliable calculation 
of K — > %% decay amplitudes in general and attempts to reproduce the the experimental value of 
e'/e and to understand the AI = 1/2 rule (the enhancement of the AI = 1/2 amplitude by a factor 
of about 22 relative to that for AI = 3/2 transitions) in particular. The non-zero experimental value 
for e'/e was historically the first evidence for direct CP- violation. In the past lattice calculations 
have tried to estimate the matrix elements by combining computations of K — > % and K — > vacuum 
matrix elements with the lowest order terms in the chiral expansion and I start with a brief review 
of the status of such calculations before proceeding to a discussion of the direct evaluation of the 
matrix elements with two-pion final states. 

4.1 K — >■ TiTi decay amplitudes from K — > it and K — > vacuum matrix elements 

At lowest order in the SU (3) chiral expansion one can determine the K^-tzti decay amplitudes 
by calculating K — > % and K — > vacuum matrix elements. In 2001, two collaborations published 
sone interesting quenched results on non-leptonic kaon decays in general and on the AI=\I2 rule 
and e'/e in particular (see table |3]). In spite of the limitations of these calculations, the authors did 
achieve the control of the ultraviolet problem, i.e. the numerical subtraction of power divergences 
and the renormalization of the weak operators. This is highly non-trivial. 
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Figure 2: The mass dependence of the K — > % matrix element. m\, m x and m z are the masses of the sea and 
valence light quark and the valence "heavy" (strange) quarks respectively. 



The RBC-UKQCD collaboration have repeated the calculation in the pion-mass range 240- 



420 MeV [43, 44]. For illustration consider the determination of (X2i, the lowest order LEC for the 
AI = 3/2 (27,1) operator: 



(27,i) = ( S( 1)l { {uu) L - {dd) L ) + (su)l (ud) L . 



(4.1) 



Satisfactory fits for the mass dependence were obtained using NLO SU(3) ChPT, but the corrections 
were found to be very large, casting serious doubt on the approach. This is illustrated in fig^| where 
the dashed curve represents the results in the SU(3) limit (m Ll = = m s ) and the value of this curve 
in the chiral limit is much below the data points, demonstrating that the one-loop corrections are 
very large. Thus the use of soft-pion theorems is not sufficiently reliable and K — > nn matrix 
elements have to be computed directly. At this conference we have also heard about a proposed 
method to combine K — > % matrix elements with ones limited to the two pions at rest [|3|]. 

4.2 Direct Calculations of K — > nn Decay Amplitudes 

From the above discussion we conclude that we must calculate K — > Tin matrix elements di- 
rectly and the RBC-UKQCD collaboration is now undertaking a major study. The ultraviolet prob- 
lem of the subtraction of power divergences (i.e. of terms which diverge as inverse powers of the 
lattice spacing) remains tractable at the expense of significant statistical errors. The infrared prob- 
lem of extracting the spectrum and amplitudes in a finite Euclidean volume is also understood as 
long as we neglect the inelastic contribution (i.e. rescattering into states other than two pions). 
Assuming that the rescattering is dominated by the s-wave states, the quantization condition for 
two-pion states in a finite volume derived by Liischer takes the form [|46|-|4"8|] 



8(q*) + <l>{q*)=nn 



(4.2) 



where 8 is the physical s-wave 7171 phase shift, is a known kinematic function and n is an integer. 
The relative momentum q* is related to the two-pion energy E by E 2 = 4(m| + q* 2 ). For illustra- 
tion, all the formulae in this section are presented in the centre-of-mass frame, but they have also 
been generalised to moving frames p^— 51]. Thus from the measured values of the energies one 
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can determine the phase-shift. The relation between the measured Euclidean matrix elements and 



the physical amplitudes is given by J52| , |53[ ] 



mi 



\A\ 2 = SnV 2 {8'(q*) + $ P \q*)} \M\ 2 , (4.3) 
*7 



where the ' denotes differentiation w.r.t. q* . (4.3) has been generalized to moving frames in [50, pip. 



4.2.1 K ->■ (7ix)[ =2 Decays 

The calculation of decay amplitudes into two pion states with isospin 2 is relatively straight- 
forward; there are no power divergences to subtract nor any disconnected diagrams to evaluate. 



We will see in section 4.2.2 that these are significant difficulties when evaluating AI = 1/2 decay 



amplitudes. An exploratory quenched study with improved Wilson fermions was completed in 



2004 [54], but at that time we had not understood the finite- volume effects at non-zero total mo- 
mentum. Results from RBC-UKQCD's exploratory quenched study with Domain Wall Fermions 
were presented in 2009 [|5| and this year M.Lightman presented results from dynamical simu- 
lations with almost physical pions, but on a course lattice [|56|]. Before discussing the results I 
mention some theoretical points. 

There are three AI = 3/2 operators whose matrix elements need to be evaluated: 

°(27,1) = { (« ; V) L - (dW) L } + (iV)i (4.4) 

3 7 /2 = (?J) L { (kV) r - (JW)*} + (jV)l {uUi) R (4.5) 
3 g /2 = (sU^l { (u j u% - + (u j d% , (4.6) 

where i and j are colour indices and (q\q2)L,R = qi7^(l T T 5 )^- The subscript (27, 1) in ([O]) 
denotes that the operator transforms as the (27,1) representation of the SU(3)l x SU(3)r chiral 
symmetry group. Oj^ are the 1 = 3/2 components of the electroweak penguin operators labelled 
C?7 8 in the standard notation for the AS = 1 effective Hamiltonian and transform as the (8,8) 
representation of SU (3)l X SU (3)r. 

A significant simplification in the calculation of the matrix elements of these operators is the 
use of the Wigner-Eckart theorem to relate the physical K + — > % + 71° matrix elements to unphysical 



K + ->■ 7T+7T+ ones [57]: 



/=2(7r + (pi)7r b2)|O^|^ + ) = -(7r + ( Pl )7r + ( K )|O^|^ + ), (4.7) 

where the subscript on the operators indicates AI Z , the change in the z-component of isospin. 
Eq. (fO) is an exact relation in the isospin limit m u = win, so that the evaluation of the matrix 

3/2 3/2 

elements of the operators Oy 2 and Oy 2 are equivalent, but there are a number of advantages in us- 

3 /2 

ing the fully extended \ k + k + ) states. The flavour structure of the operators Oy 2 is simply (sd){ud) 

rather than the (sd) iuu — dd) + (su)ud structure of the operators Oy 2 m ec l- Q4.4| ) - d4.6| ). 

With the 7T + 7r + -state we can impose antiperiodic boundary conditions on the J-quark say, so 
that the ground state is \n + (n / L)n + (—n / L)) , where the arguments represent the momenta in lat- 
tice units (up to finite-volume corrections). In contrast to the use of periodic boundary conditions, 
for which the ground state corresponds to the two-pions at rest, it is not now necessary to isolate an 
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Figure 3: Schematic diagram illustrating the contractions for AI = 3/2 K — > nn decays in the RBC-UKQCD 
calculation [|55|]. 



excited state in order to have a decay into two-moving pions. For the physical decay, the minimum 
size of the lattice is halved from about 6 fm to 3 fm. 

The final theoretical point I wish to make here is that the use of the Wigner-Eckart theorem 
also allows us evaluate the Lellouch-Luscher factor relating the measured matrix elements to the 



physical amplitudes directly [p8||. In particular, this factor requires knowledge of the derivative 



of the phase-shift (see eq. (|4.3|)). By imposing partially-twisted boundary conditions [|59J-|62J on 
the <f-quark with twisting angle 6 (it is sufficient to perform the twist in a single direction), the 
two-pion ground state now corresponds to a pion with momentum d/L and the second pion with 
momentum (6 —2n)/L. The corresponding %% s-wave phase-shift can then be obtained by the 
Liischer formula Q4.2| ) as a function of 6 which allows for the derivative of the phase-shift to be 
evaluated directly at the masses being simulated. This procedure was tested in an exploratory 
simulation [ |58| ] but has not yet been implemented in the main RBC-UKQCD programme. 

I now briefly summarise the results from the RBC-UKQCD Collaboration presented at this 
conference [56]. The correlation functions which need to be evaluated are illustrated in fig.||. The 



simulations were performed on 32 3 x 64 x 32 lattice using a Domain Wall Fermion action for the 
quarks and the DSDR (dislocation suppressing determinant ratio) gauge action on a course lattice 
(cr 1 ~ 1.4 GeV). The motivation for such a coarse lattice was to enable almost physical pions to 
be simulated, the K — » %% amplitudes were obtained with a partially quenched pion of mass m n = 
145(6) MeV (the corresponding unitary pion has m n ~ 180MeV) and the kaon mass is 519(2) MeV, 
both close to their physical values (m % = 139.6 MeV, = 493.7 MeV). The properties of these 
ensembles were discussed at this conference by Bob Mawhinney 

With the masses given above {m % = 145.6(5) MeV, = 519(2) MeV) and choosing the mo- 
mentum of the pions to be y/2n/L, the energy of the two-pion system is found to be 516(9) MeV so 
that the kinematics is almost matched. Sample plateaus for each of the three operators are presented 
in fig. The preliminary result for the real part of the amplitude is ReA2 = 1.56(07) stat (25) sys t x 
10 GeV, to be compared to the experimental number of 1 .50 x 10~ 8 GeV . The non-perturbative 
renormalization has not been completed for the electroweak operators which contribute to Im(A2); 
guesstimating the renormalization factors RBC-UKQCD quote Im(A 2 ) = — (9.6 ± 0.04 ± 2.4) x 
10~ 13 GeV [56], but after the renormalization is complete the systematic error will decrease from 
about 25% to 15%, dominated by the uncertainty in the value of the lattice spacing. It appears that 
these calculations are possible with good precision. 

4.2.2 K -)■ (7r7r)/ = o Decays 

The evaluation of AI = 1/2 matrix elements is very considerably harder than that for AI = 
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Figure 4: Sample plateau for the matrix elements (a) {scI)l (ucI)l, (b) (sd)i (ucI)r and (c) (s'd-')i (uJd l )R, 
where (sd)i (M)l.r — {s'y^(l — J 5 )d') («- 7 7 M (l T Y 5 )d-'), and i,j are colour labels. 




Figure 5: The four diagrams which contribute to the two-pion propagator. 



3/2 operators. The two-pion state with 1 = has vacuum quantum numbers and the vacuum 
contributions have to be subtracted requiring large statistical cancelations to obtain the exp(— E n% t) 
behaviour, where E nn is the energy of the two-pion state. To illustrate this, consider the four 
diagrams in fig|5] which make up the two-pion propagator. The 1 = 2 %% correlation function is 
proportional to D-C whereas the 1 = correlation function is proportional to 2D+C-6R+3V. The 
major practical difficulty is to subtract the vacuum contribution with sufficient precision. 

To demonstrate the above, I now present some results from the exploratory study by RBC- 
UKQCD on a 16 3 x 32 x 16 lattice with an inverse lattice spacing of cr 1 = 1.73 GeV and an 
unphysically heavy pion with mass 420 MeV [|64|]. To increase the statistical precision, quark 
propagators are generated from each of the 32 time slices. The four components of the 1 = two- 
pion correlation function 2D+C-6R+3V are shown separately in fig.[6| []64|], from which we see that 
the error on the 3V component grows significantly at larger times. 

In spite of the difficulty of subtracting the vacuum contributions, at this conference Qi Liu 
presented the results of a complete calculation of Aq and A 2 , albeit at unphysical kinematics with 
m K ~ 420 MeV and with the pions at rest [54]. There are 6 diagrams, shown in fig.|7| and a total 
of 48 Wick contractions. The diagrams labeled by Mix3 and Mix4 in fig.[7| are needed to subtract 
the power divergences which are proportional to matrix elements of the pseudoscalar density sy=,d. 
The results for this unphysical kinematics are [|4j 



Re A = (3.0 ±0.9) 10~ 7 GeV, 
Re A 2 = (5.394 ± 0.045) 10~ 8 GeV, 



Im A 
Im A 2 



-(2.9±2.2)10~ n GeV, 



(4.8) 



(7.79 ±0.08) 10 GeV. (4.9) 



The precision is limited by the disconnected diagrams and the vacuum subtraction. At least the 
calculation of the real part of Ao appears to be tractable for this simplified kinematics at heavy 
pion masses, whereas for ImAo we will require more statistics to confirm that there is indeed a 
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Figure 6: Contribution of the four components to the 7 = two-pion correlation function. 
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Figure 7: Types of correlation function corresponding to K — > Tin decays. The blue circles represent the 
insertion of a four quark operator appearing in the weak Hamiltonian and the squares that of the pseudoscalar 
density. 



signal. The greater challenge is now to proceed towards performing the corresponding calculations 
at physical kinematics. 



5. r\ and r\' Mesons 



In the previous section we saw the importance of controlling disconnected diagrams and this is 
the case for many important phenomenological quantities. Here I mention one other longstanding 
issue, the spectrum and mixing of tj and v\' mesons. To determine the propagators we need to 
evaluate the diagrams shown in fig.|8| and I report here on our recent study on a lattice of spatial 



extent 16 and lattice spacing a =1.73 GeV [65]. Using interpolating operators of the form 



Oi = — — ' and O s = sy 5 s, 
v2 



(5.1) 
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Figure 9: (a) Values of each of diagrams as a function of t. (b) The effective mass plots for the rj and rj' 



mesons. In each case the results are for the lightest mass in the simulations; see 1 65 ] for details. 



we calculate the correlation functions 

1 31 

M0 = ™ £ ( O a (t+t')O p (t') ) where a,P=l,s (5.2) 
51 t'=0 

= (Cu-2D n -V2D ls \ 
\-y/2D a C ss -D ss ) 

where the diagrams are defined in fig.|8| The propagators are generated with sources at each of 
the 32 time slices t' . In fig.|9](a) we present the values of the diagrams at the lightest mass used in 
the simulations as a function of t, from which we see that the usual expectation that disconnected 
diagrams are small does not apply here and also the expected feature that the uncertainties on the 
disconnected diagrams are significantly larger than those on the connected ones. The effective 
mass plots for the eigenvalues of X(t) are shown in fig. §(b), again at the lightest masses used in the 
simulations, from which we see a good plateau for the tj and a short, but clear, plateau for the tj'. 
After the chiral extrapolation we obtain the values = 573(6) MeV and = 947(142) MeV, 
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where only the statistical errors are shown, compared to the experimental values of 548 MeV and 
958 MeV respectively. 



As explained in ref. [65], it is also possible to determine the mixing angle. Once SU(3) flavour 
symmetry is broken it is not clear that the physical v\ and v\' states are simply linear combinations 
of the octet and singlet flavour combinations; nevertheless it is a standard phenomenological as- 
sumption that they are. Based on this assumption we obtain a value of 6 = —14.1(2.8)° (statistical 
error only) in agreement with the phenomenological estimates lying in the range —20° to —10°. A 
particularly important feature of the calculation is that it is possible to check the orthogonality of 
the mixing matrix confirming that it is a good approximation to consider only the symmetric octet 
and singlet states in order to understand the mixing. 

In spite of the limited precision, our calculation demonstrates that QCD can explain the large 
mass of the ninth pseudoscalar meson and its small mixing with the SU(3) octet state and provides a 
first benchmark for future calculations. There remains much to be done now to reduce the statistical 
errors and to quantify in detail the systematic uncertainties. Progress in the efficient determination 
of all-to-all propagators will be very important in improving the reach and precision of lattice 
calculations of physical quantities in which disconnected diagrams play an major role. 

6. Long-Distance Contributions to Physical Quantities 

We are used to calculating the short-distance contributions to physical processes, generally 
formulating the calculation as the evaluation of the matrix element of a local operator. A good 
example is the evaluation of the Bk parameter in neutral-kaon mixing illustrated by 

s -> ? d s 




W? >W = C(Miy/jU) 

d s d 

The determination of the matrix element of the resulting AS = 2 local operator has been discussed 
in the section ||. In many cases the short-distance contribution is the dominant term, but long- 
distance contributions are not always negligible, for example if the GIM suppression is logarithmic 
or if there is a CKM enhancement (even if the GIM suppression is power like). As lattice results 
become more precise, we should try to compute the long-distance contributions effectively and this 
represents a new type of calculation. Early thoughts in this direction include ref. for rare kaon 
decays and ref. [|57|] for neutral kaon mixing. 

6.1 Rare Kaon Decays 



In ref. [|66|], the authors conclude that it is possible in principle to evaluate the long distance 
effects for K — > 7il + £~ and K — > nvv decays. This requires the evaluation of r -products of the 
form 

T Q .j(q 2 ) =N V f d A xd A ye-^ (n\T [Q(x)J^y)} \K) (6.1) 

where Q is one of the four-quark operators appearing in the weak Hamiltonian, J is a weak or 
electromagnetic current and Ny is a volume factor. The generic lattice calculation is therefore of 
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correlation functions of the form 

-ij d 4 xe-^ x {0\<l> n (^,p)4{x) [fi"(0)-fif(0)] 4(t K ,k)\0), (6.2) 

with tn > and t K < 0. § n and 0^- are the interpolating operators for the % and K mesons respec- 
tively (after the three-dimensional Fourier transform has been taken). The main issue discussed in 
ref. [^] is that of renormalization, the subtraction of power divergences and the consequences of 
contact terms. The authors conclude that, as a result of symmetries and the GIM mechanism, the 
power divergences can be removed and they check their general arguments by one-loop perturba- 
tive calculations. They believe that their study opens a new field of interesting physical applications 
for the lattice community, although to date no such numerical calculations have been performed. 

6.2 Long-distance contribution to the Kl-Ks mass difference and Ek 

At this conference Norman Christ presented some interesting ideas how the long-distance 



contribution to the Kl — K$ mass difference and £k might be evaluated [67]. This requires 



(i) the calculation of the long-distance contribution to the matrix elements of the product of two 
AS = 1, four-quark weak operators between kaon states; 

(ii) the subtraction of the short distance part of this matrix element in a way that is consistent 
with the original explicit evaluation of the short-distance contribution; 

(iii) a generalization of the Lellouch-Luscher approach to finite-volume corrections to second 
order in the weak interaction. 

Although much work remains to be done to develop these ideas into a practicable method, the main 
theoretical steps have now been taken. 



7. Non-leptonic 5-Decays 

I end this talk with a discussion of an important class of processes in flavour physics for which, 
up to now at least, little or no progress has been made in formulating a lattice approach, namely 
non-leptonic B-decays. A huge amount of precise information, from over 100 channels, has been 
obtained about decay rates and CP-asymmetries for the exclusive decays of B-mesons into two light 
mesons. Unfortunately, with just a few exceptions (most notably the CP-asymmetry in the golden- 
mode B — > JyKs), our ability to deduce fundamental information about CKM matrix elements is 
limited by our inability to quantify the non-perturbative QCD effects sufficiently precisely. 

The operator product expansion can be used to separate the long and short distance physics, 
so that the non-perturbative QCD effects are contained in the matrix elements of local operators, 
(MiM2|O;(0)|S), where 0, is a AB = 1 operator and are light mesons. In contrast to K — > 
7171 decays discussed in section |], where it is a good approximation to restrict consideration of 
final-state rescattering to two-pion states, the large mass of the B-meson means that very many 
intermediate states contribute. We do not yet(?) have a theoretical framework for treating this in 
simulations in Euclidean space. 
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Figure 10: Illustration of nave factorization for the decay B ( [ — > n + 7i~ . The local four-quark operator 
(ub)v-A (du)v-A is shown as a product of the bilinear operators (ub)v-A an d {du)v-A represented by the 
black circles. The two black circles are separated for clarity. 



In the past phenomenological approaches to nonleptonic B-decays were based on naive fac- 
torization in which the matrix element was reduced to products of matrix elements of B — > M\ and 
vacuum—^ M2 and/or vice- versa, depending on the flavour quantum numbers. This is illustrated in 
fig.[K| where the matrix element of the local four-quark operator (iib)v-A {du)v-A 1S shown as the 
product of the matrix elements of the two currents: 

( it + it~ I (ub) V -A {du) V - A \B d ) = (d~u) V - A 1 ) ( 7T+ 1 {ab)v- A \B d ). (7.1) 



The motivation for assuming (7.1) is that the two factors on the right-hand side are known or 
calculable (in principle at least). The first factor is proportional to the leptonic decay constant f K 
and the second is proportional to the form-factors of semileptonic B — > % decays. On the other hand 
the limitations of such a model are clear; the renormalization-scale dependence of the two sides do 
not match and rescattering effects are not included. 

In 1999, together with Beneke, Buchalla and Neubert, we realised that as my — > °°, the long- 



distance QCD effects do factorise into simpler universal quantities Q68|— [70(] 

{M u M 2 \Oi\B) = £/f^ Ml (m2) J^duT l I j {u)<S> M M + (Mi M 2 ) 



j 



+ / rf£dudvr/ 7 (£,a,v)0i l (£)0M 1 (v)0 i fc(w), (7.2) 
Jo 

where fJ^ Mi are the form factors for B — » M\ transitions, the <J> are light-cone distribution ampli- 
tudes and T I,n are short distance contributions and are calculable in perturbation theory, u and v are 
the momentum fractions carried by the quarks in the mesons. The significance of the factorization 
formula stems from the fact that the non-perturbative quantities which appear on the right-hand 
side are much simpler than the original matrix elements. They either reflect universal properties 
of a single meson state (the light-cone distribution amplitudes) or refer to a B — > meson transition 
matrix element of a local current (form factors). Conventional (naive) factorization is recovered as 
a rigorous prediction in the infinite mass limit (i.e. neglecting 0(a s ) and 0(A QCD /mb) corrections). 
Perturbative corrections to naive factorization can be computed systematically and the results are, 
in general, process dependent. It is a remarkable feature that all strong interaction phases are gener- 
ated perturbatively in the heavy quark limit. The factorization formulae are valid up to 0{K QCO /mb) 
corrections and the main limitation of the framework is due to the fact that since m\, is not so large, 
CKM and chiral enhancements to non-factorizable 0(AQCD/mb) terms are important. 
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Although we do not know how to evaluate the matrix elements for B — > M1M2 decays directly, 
we can ask what can lattice simulations contribute to the factorization formula ( [7.2| )? The moments 
of the light-cone distribution functions for the light mesons can an d are being computed [fTl]] 
as are the B — > M form factors. What we do not know how to compute at this stage are the parton 
distribution amplitudes of <J># or its moments and I end this section with a brief explanation of the 
reasons for this. <J># is defined by 

*Bap(k + ) = J dz-e^-(0\up( Z )[z,0}ba(0)\B)\ z+zx=0 , (7.3) 

where ± denote light-cone coordinate and [z,0] represents the part-ordered exponential of gauge 
fields between z and 0. In evaluating matrix elements, <J># is convoluted with the perturbative 
hard-scattering amplitude T t n and the relevant quantity is 

y/2 r dk+ ,~ . 

/ T ±0 fl (A+). (7.4) 
H Jo k + 

(In higher orders of perturbation theory factors containing \og(k + ) appear.) Although at large k + , 
<fe(£+) ~ V^+> the convolution in eq.(|L4|) is finite. In lattice calculations, at least up to now, we 
know how to calculate the matrix elements of local operators and the positive moments of 0# (k + ) 
can indeed be written in terms of local operators. However they diverge as powers of I /a and we 
still need to develop techniques to subtract these divergences with sufficient precision. 

This discussion underlines the fact that we need new theoretical ideas for lattice simulations to 
start contributing to the evaluation of the non-perturbative QCD effects in B — > M1M2 decays and 
hence to enable fundamental information to be obtained from the wealth of experimental measure- 
ments of the rates and CP-asymmetries. 



8. Conclusions 



At this conference we have seen many beautiful contributions to particle physics phenomenol- 
ogy, both in improved precision and in the extensions of computations beyond the standard quan- 
tities. We readily forget that it was only a few years ago that results presented at the annual lattice 
symposium were largely in the quenched approximation, with an error which was not possible 
to quantify reliably. We then moved on to a brief period with dynamical quarks with masses of 
0(500) MeV or so until today we arrive at simulations with almost physical pions. This improve- 
ment has to be continued vigorously if precision flavour physics, which has been the focus of this 
talk, is to play a complementary role to large p± discovery experiments at the LHC in unravelling 
the next level of fundamental physics. As was noted in sec. [7] however, there are quantities for 
which a large amount of experimental data is available, and yet for which we do not yet know how 
to begin formulating the calculations to make them accessible to lattice studies. 

At the previous lattice conference which Guido Martinelli helped to organise which was held 
in 1989 in Capri, Ken Wilson made the seemingly pessimistic prediction that it will take about 30 
years to have precision Lattice QCD. We only have 9 years left to fulfill the prediction, but we are 
now well on our way. 
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